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ABSTRACT 

We report the identification of a very metal-poor damped Lya system (DLA) at z a b s = 
3.067295 that is modestly carbon-enhanced, with an iron abundance of ~ 1/700 solar 
([Fe/H] = -2.84) and [C,0/Fe] +0.6. Such an abundance pattern is likely to be the result of 
nucleosynthesis by massive stars. On the basis of 17 metal absorption lines, we derive a 2cr 
upper limit on the DLA's kinetic temperature of Tola < 4700 K, which is broadly consistent 
with the range of spin temperature estimates for DLAs at this redshift and metallicity. While 
the best-fitting abundance pattern shows the expected hallmarks of Population III nucleosyn- 
thesis, models of high-mass Population II stars can match the abundance pattern almost as 
well. We discuss current limitations in distinguishing between these two scenarios and the 
marked improvement in identifying the remnants of Population III stars expected from the 
forthcoming generation of 30-metre class telescopes. 
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1 INTRODUCTION 

With the advent of the Sloan Digital Sky Survey in combina- 
tion with 8-10 metre class telescopes, it has become possible 
to routinely identify clouds of gas in the high redshift Universe 
that h arbour the metals from some of the earliest stellar popula- 
tions jPettini et al]2008al : IPenprase et al]201fi iBecker et al]201 ll ; 
ICookeetal .11201 lbh . The challenge we are now faced with is to 
measure the allowable ranges of mass, explosion energy, and metal- 
licity of these early stars, using metal-line diagnostics alone. In the 
last decade, attention has been drawn to the diagnostic power of 
carbon for probing the nature of zero or near-zero metallicity stars. 
In particular, it is believed that metal-free stars are more likely 
to yield an abundance pattern th at is enhanced in carbon relative 
to iron dUmeda & Nomotdl2003h . In fact, such signatures may be 
observed among the population of metal-poor stars in the halo of 
our Galaxy that exhibit a marked carbon-enhancement in their at- 
mosphere; these stars are known as carbon-enhanced metal-poor 
(CEMP) stars, and are currently defined to have [C/Fe] > +0.7 



* Based on observations collected at the W.M. Keck Observatory, which is 
operated as a scientific partnership among the California Institute of Tech- 
nology, the University of California and the National Aeronautics and Space 
Administration. The Observatory was made possible by the generous finan- 
cial support of the W. M. Keck Foundation, 
t email: rcooke@ast.cam.ac.uk 



dAoki et al] 12007. although some authors adopt [C/Fe] > +0.5 or 
> +1.0 as the defining cut)[] 

Several possibilities have been put forward to explain the pe- 
culiar abundance patterns of these stars. One option is that they re- 
sult from mass transfer from a now extinct asymptotic giant branch 
companion star. This type of CEMP star also exhibits enhance- 
ments in its s-process neutron-capture elements, and is thus labelled 
a CEMP-s star. The above mecha nism is now a well-e stablished 
means to produce a CEMP-s star dLucatello et~ai1l2005h . Alterna- 
tively, the carbon enhancement may be the residual signature from 
a previous generation of stars — possibly Population III — that 
seeded (w ith high carbon ab undance) the birth cloud of the star we 
see today dRvan et al Jl 20051) . These are known as CEMP-no stars, 
as they show no strong enhancement in their neutron-capture ele- 
ments. Perhaps the strongest empirical evidence in support of the 
latter picture is the observation of an increasing fraction of C EMP 
stars with decreasing metallicity l lBeers & Christliebll2005l) . Fur- 
thermore, the fraction of CEMP stars that are labelled as CEMP-no 
stars becomes rel atively more common at the lowest metallicities 
dAoki et alj|2007l) . 

Additional evidence in support of this picture has recently 
been provided by the discovery of a CEMP damped Lya- system 
(DLA) at z abs = 2.340 with [Fe/H] * -3 and [C/Fe] = +1.53 

1 We adopt the standard notation: [A/B] = log(Ar A /W B ) - log(iV A /iVB)o, 
where Na.b refers to the number of atoms of element A and B, and the 
second term refers to the solar ratio. 
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l lCookeetal]|2011<J) . These authors proposed that such DLAs 
could be the missing link between the enrichment of a primor- 
dial cloud by Population III stars, and the later incorporation of 
the trace metals into objects such as the CEMP-no stars seen in the 
halo of our Galaxy. Indeed, models of metal-free nucleosynthesis 
are able to suc cessfully reproduce the observed abunda nce pattern 
for this DLA l lKobavashi. Tominaga. & Nomotolfeoill) . However, 
there may be an additional contribution to the metal content of this 
DLA from long-lived second ge neration asymptotic giant branch 
stars dSalvadori & Ferrarallioij) , which are known to be prolific 
carbon producers. 

DLAs are clouds of neutral gas seen in absorption against the 
light of backgroun d QSOs; their prop erties h a ve bee n reviewed by 
IWolfe et all d2005h . As dis cussed by IPettinil fc004h and more re- 
centlv lCooke et alj feOlibh . measurements of their chemical com- 
position complement very effectively local stellar abundance stud- 
ies. While most DL As are metal-poor, with typical metallicities 
of ~ 1/30 of solar fellison et al]|2012h . interest has recently fo- 
cused on identifying and studying the most metal-poor DLAs, 
with [Fe/H] < -2, which may still retain the nucleosynthetic sig- 
nature s of the first episodes of metal en r ichment in the Uni- 
verse jErni etalJ 120061 ; IPettini et alj l2008al : IPenprase etai]|2010i ; 
ICookeetalfeOllbl) . 

Empirically, it is found that the most metal-poor DLAs tend 
to have very quiescent kinematics, with most of the metal line 
absorption concentrated in one or two components of low ve- 
locity dispe rsion. For example, the z a bs = 2.340 CEMP DLA 
reported by ICooke etaf consists of a single absorption 

component with an internal one-dimensional velocity dispersion 
cr = 1.7 km s -1 . In this extreme regime it becomes important to 
differentiate between bulk gas motions which broaden all absorp- 
tion lines by the same amount, and thermal motions which, for the 
same gas temperature, broaden lighter elements more than heavier 
ones (see Section[3). Unfortunately, with present instrumentation it 
is generally not possible to fully resolve the metal line profiles pro- 
duced by such quiescent gas, although under favourable circum- 
stances a curve-of-growth analysis can help to separate the rela- 
tive contributions of t urbulent and thermal broadening to the line 
equivalent widths (e.) 



iTumlinson et al. 20K 



O'Meara et al 



Carswell et al 



200ll ; |jorgenson et alj2009l ; 



201 lh . 



In their analy sis of the CEMP DLA at z abs = 2.340 
ICooke et al.U2011af) assumed that turbulent motions dominate, as is 
the case if the gas temperature is less than a few thousand degrees. 
Higher temperatures would have the effect of decreasing selectively 
the carbon abundance, and hence the degree of carbon enhance- 
ment (see Section[3]for further detai ls), although the resulting line 
profil es are a poorer fit to the data dCookeluOl id ; ICarswell et alj 
l2012h . 

In order to make progress on these issues it is important to 
identify and study further examples of CEMP DLAs. In this paper 
we report the discovery of a new candidate with a lower carbon en- 
hancement than the first example of this class, but where the gas 
temperature can be constrained within a narrower range. In Sec- 
tion|2]we describe the observations and data reduction procedures, 
saving the discussion of line profile fitting and the abundance anal- 
ysis to Section[3] We discuss the possibility of using such systems 
to measure the kinetic temperature of DLAs in Section [4] before 
comparing the best-fitting abundance pattern to models of stellar 
nucleosynthesis in Section[5]and concluding in Section[6] 



2 OBSERVATIONS AND DATA REDUCTION 

To search for additional candidate CEMP DLAs out of the increas- 
ing database of damped systems, we target metal-poor systems 
showing an unusually strong Cii/11334 absorption line in spec- 
tra recorded at medium resolution (full width at half maximum, 
FWHM=: 50-100 km s~'). While this relatively coa rse resolving 
power is insufficient for a reliable abundance analysis dCooke et alj 
1201 lbl) . it does provide the means to identify some of the most 
promising candidates. 

We selected the DLA at z a b s = 3.0675 towards the quasar 
J1358+6522 on the basis of the Echelle Spectrogr aph and Imager 
(ESI) observations with FWHM = 60 km s~' by IPenprase et alj 
l 2010h . These authors reported [O/H] = -3.08 ± 0.15 and [C/O] = 
+0.44±0.43; given that DLAs in this metallic ity regime commonl y 
ex hibit an a-enhancemen t of [O/Fe] - +0.4 dCooke et al]|2011bl) . 
the lPenprase et alj J20T0I) data suggest that [C/Fe] ^ +0.8. 

We conducted follow-up observations of J1358+6522 for 
28 950 s divided into nine exposures with the W. M. K eck Obser- 
vatory 's High Resolution Echelle Spectrometer (HIRES; IVogt et all 
1 19941) on 2011 March 23, 24 (program ID: A152Hb). We used the 
C5 decker (a 7.0 x 1.148 arcsec slit) which delivers a nominal res- 
olution of 8.1 kms~' for a uniformly illuminated slit (but see Sec- 
tion[3] where we adopt a lower FWHM resolution of 7.0 km s due 
to the sub-arcsecond seeing conditions), which we sample with ~ 3 
pixels. We used the red cross-disperser with the WG360 filter and 
2x2 on-chip binning. 

The data were reduced with the makee data reduction pipeline 
maintained by Tom Barlow)^ makee performs the usual steps rel- 
evant to echelle data reduction, including bias subtraction, flat- 
fielding, order definition and extraction. Each science exposure was 
followed by a ThAr hollow-cathode lamp frame which was used 
to wavelength calibrate the extracted orders. The wavelength scale 
was converted to vacuum heliocentric using the software package 
uves_popler 0This software was also used to combine the extracted 
orders from all exposures, whilst rejecting deviant pixels and un- 
usable orders. Finally, the combined spectrum was normalised by 
fitting and dividing out the quasar continuum and emission lines. 
The metal absorption features of interest were then extracted in 
±200 km s" 1 windows about the line centroid, and a fine adjustment 
to the continuum was applied if necessary. The resulting signal-to- 
noise ratio of the data at 5000 A is S/N ~ 30 per resolution ele- 
ment. Examples of absorption lines in the z abs = 3.0675 DLA are 
collected in Figure Q] 



3 ANALYSIS AND RESULTS 
3.1 Profile Fitting 

As can be seen from Figure Q] and Table Q] our HIRES spectrum 
covers many transitions of the dominant ionisation stages of C, N, 
O, Al, Si, S and Fe, as well as lines in the Lyman series of Hi. 
All the metal lines exhibit a single, narrow absorption component 
which, with FWHM = 5 pixels, is just resolved with our instrumen- 
tal resolution of ~ 3 pixels. 

In order to deduce the corresponding element abundances, 
we modelled the absorption lines with theoretical Voigt profiles 



2 We used makee version 5.2.4, available for download from 
http : //spider . ipac . caltech . edu/staf f /tab/makee/ 

uves_popler can be downloaded from 
http : //astronomy . swin . edu . au/ -mmurphy/UVES.popler/ 
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Figure 1. Selected absorption lines in the HIRES spectrum of the Zabs = 3.067295 DLA towards the QSO J1358+6522. Black histograms are the data, and 
red continuous lines show the best fitting theoretical profiles computed with vpfit as discussed in Section [5] The top panel shows the Hi Lyo- absorption line 
together the damped profile corresponding to a column density log |W(H i)/cm~ 2 ] = 20.47 ± 0.07. Blue lines in some of the panels are used to indicate unrelated 
Lya forest blends. The weak feature at +26 km s in the Si n /il260.4221 panel is Fe n ,11260.533. The y-scale is residual intensity. 



(convolved with the instrumental broadening function) using the 
X 2 -minimisation software vPFrfl which returns the best-fitting 



Doppler parameter b = V2rr and redshift z a t, s of the DLA, in addi- 
tion to the column densities N of the available ions. The line fitting 
procedure consisted of three main steps, which we now outline. 



4 vpfit is available from 

http : //www . ast . cam . ac . uk/~rf c/vpf it .html 
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Table 1. Ion column densities of the DLA towards J1358+6522 at z aDS = 3.067295 



Ion 


Transitions used 


log N(X)/cm 2 


l0 § £ ( X >DLA 


log e(X)°' c 


fX/H]*-' 7 


[X/Fe] 


Hi 


1025, 1215 


20.47 ± 0.07 




12.0 






Cii 


1036, 1334 


14.65 ±0.10 


6.18±0.10 


8.43 


-2.25 ±0.10 


+0.59 ±0.10 


Ni 


1200.2 


12.62 ±0.14 


4.15 ±0.14 


7.83 


-3.68 ±0.14 


-0.84 ±0.14 


Oi 


971, 976, 988x3, 1039, 1302 


14.94 ± 0.05 


6.47 ± 0.05 


8.69 


-2.22 ± 0.05 


+0.62 ± 0.06 


Aln 


1670 


11.92 ±0.03 


3.45 ± 0.03 


6.44 


-2.99 ± 0.03 


-0.15 ±0.04 


Sin 


989, 1260, 1304, 1526 


13.40 ± 0.03 


4.93 ± 0.03 


7.51 


-2.58 ±0.03 


+0.26 ± 0.04 


Sn 


1253, 1259 


13.11 ±0.09 


4.64 ± 0.09 


7.14 


-2.50 ± 0.09 


+0.34 ± 0.09 


Fen 


1144, 1260, 1608 


13.10 ±0.03 


4.63 ± 0.03 


7.47 


-2.84 ± 0.03 


+0.00 



"log e(X) = 12 + log JV(X)//V(H). 

b We have not included the uncertainty in H i. 

lAsplund et a?] 120091) . 

d [X/H] DLA = loge(X) DLA - loge(X) G . 



3.2 Instrumental Resolution 

In the first step, we determined the true value of the instrumental 
resolution which is likely to be lower than the nominal value of 
FWHM= 8.3 km s appropriate to the uniformly illuminated slit 
provided by the C5 HIRES decker, given that the seeing during our 
observations was better than the 1.15 arcsec width of the slit. Exam- 
ination of the emission line profiles from the Th-Ar hollow-cathode 
lamp showed that, at the signal-to-noise ratio of our data, a gaus- 
sian profile is a good approximation to the instrumental broadening 
function. 

We assessed the resolution of our spectrum empirically from 
the metal lines in the za, s = 3.0675 DLA — which are the nar- 
rowest features available in the spectrum — using vpfit in a se- 
ries of model fits which differed by the assumed value of the 
instrumental broadening F WHM and noting th e corresponding 
value of x 1 , as described in ICooke et al. I feOllal) . The minimum 
value of x 1 was found to be that for an instrumental resolution 
FWHM illstr = 7.0kms~', which we adopted in all the subsequent 
model fitting. Lowering the instrumental resolution from the nom- 
inal FWHM, nstr = 8.3 km s _1 , to FWHM instr = 7.0 km s" 1 has only 
a modest effect on the derived element abundances, since for most 
elements we cover at least one absorption line which is optically 
thin. Thus, for example, the change in the column density of Fe n 
is only AlogA'(Feii) = -0.02. The largest change is for Cn, 
AlogAf(Cn) = -0.07, because both Cn/11334 and ,11036 transi- 
tions are moderately saturated (see FigureQ}. 



3.3 Abundance Analysis 

Having determined the most likely value of the instrument resolu- 
tion, we used vpfit to derive the absorption parameters, or 'cloud 
model', that best represent the data. When transitions from ele- 
ments of different atomic mass are available, as is the case here, 
vpfit can solve separately for the macroscopic, or turbulent, broad- 
ening parameter, b tmb , which is independent of atomic mass, and 
the microscopic, or thermal, broadening parameter, b xb , which is 
related to the temperature T of the gas and the mass of the ion m 
via the standard relation: bi = 2kT/m, where k is the Boltzmann 
constant. The two contributions add in quadrature to give the total 
broadening: b 2 = b 2 ,+b 2 ,. 

o tot turb th 

Given the simplicity of the kinematics of the DLA under con- 
sideration, where all the absorption lines due to ions that are dom- 
inant in neutral gas (in our case Oi, Ni, Cn, Aln, Sin, Sn and 



Fen) appear to consist of a single, narrow component, we forced 
vpfit to find the best-fitting values of z aDS , feturb and T common 
to all the above species. We found that a single component with 
redshift z abs = 3.067295 ± 0.000001, Doppler parameter b tot = 
3.00 ± 0.06 km s , and the ion column densities (and associated 
errors) as listed in Table Q] provides the best representation of the 
data. The corresponding line-profile fits are overplotted on the data 
in Figure Q] Interestingly, vpfit favours low values of temperature 
such that the broadening is dominated by turbulent motions (further 
details in Section \3A[ . 

Since the species observed are the dominant ion stages of 
the corresponding elements in Hi regions, the element abun- 
dances can be obtained directly by dividing the ion column den- 
sities in the third column of Table Q] by the H i column density 
logA f (H i) = 20.47 ± . 07. C omparison with the solar abundance 
scale of lAsplund et"aH d2009h finally gives the abundance pattern 
relative to solar listed in the penultimate column of Table [T] 

We find [Fe/H] = -2.84 ± 0.03, making this DLA one of the 
most metal-poor known. We also find [O/H] = -2.22 ± 0.05 and 
[C/O] ^ ; these values are sig nificantly different from those re- 
ported by IPenprase et al .1 d2010h from data of lower spectral res- 
ol ution, [Q/H] = -3.08 ± 0.15 and [C/O] +0.45. As discussed 
bv lCookeetal] ( l2011bl) . such differences are not unexpected when 
comparing medium- and high-resolution data (and are indeed the 
reason why high resolution follow-up observations of candidate 
metal-poor DLAs are essential for accurate abundance determi- 
nations, despite the substantial investment in observing time re- 
quired). 

In Figure|2]we show the abundance pattern of the DLA relative 
to Fe, together with that of a typical 'carbon-normal' DLA, taken 
from th e recent survey of very metal-poor DLAs by ICooke et al .1 
(2011b). Both C and O appear enhanced relative to Fe by larger 
factors than is normally the case, although the C enhancement is 
much less pronounced than in the first reported case of a CEMP 
DLA, [C/Fe] = +1.53 in the z abs = 2.340 DLA in line to the QSO 
J0035-0918 dCooke et al]|201 lal) . Before considering possible in- 
terpretations of the chemical composition we have uncovered (Sec- 
tion[5), we discuss the effect of thermal broadening on the derived 
element abundances, estimated through a third and final set of vpfit 
simulations. 
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Figure 2. The abundance pattern for the carbon-enhanced DLA (black 
boxes, where the height of each box represents the l<x statistical uncer- 
tainty) is shown for the best-fitting cloud model. For comparison, we also 
show the abundance pattern of a typical very metal-poor carbon-normal 
DLA (gree n boxes, the height of each box shows the lcr dispersion in the 
population; ICooke et all201 lbl) . The dashed line is drawn at the solar value. 



3.4 Thermal Broadening? 

As explained in the previous section, the cloud model which best 
fits the observed metal line profiles is one in which thermal motions 
do not make a significant contribution to the line broadening, i.e. 
the gas in this DLA appears to be 'cold'. We deduced an upper limit 
to the gas temperature by no longer treating b th as a free parame- 
ter in the line fitting procedure, but instead forcing vpfit to adopt a 
fixed value of T (foturb remained a free parameter th roughout these 
tests). In a pr ocedure similar to that described by ICarswell et al.l 
j2012l see also lLampton. Margon. & BowveJ 19761) . we considered 
a grid of temperatures within the range 100 K < T < 1.3 x 10 4 K, 
and for each value of T, rerun vpfit to find the best-fitting model 
parameters and corresponding value of x 2 ■ m these tests we ex- 
cluded the weak N 1 and S 11 lines (see Figure[TJ, because they offer 
no constraint on the line broadening. 

In Figure [3] we plot the resulting differential x 2 value, 
A^C^dla) = x 2 ( t dla) - xLn for each temperature, where x 2 mn 
refers to the x 1 of the best-fitting model. Noting that seven param- 
eters can alter the total line width (i.e. the thermal and turbulent 
Doppler widths in addition to five column density estimates which 
are weakly anticorrelated with the total line width), the 70 and 95 
percent confidence regions of the x 2 distribution (i.e. lcr and 2cr 
for a normal distribution) correspond to Ax 2 = 8.4 and 14.1 re- 
spectively. These are indicated with dashed lines in Figure[3] From 
this test we derive a 2cr upper limit to the DLA's kinetic tempera- 
ture of Tola < 4700 K. Furthermore, we remark that a cloud model 
where the line profiles are broadened entirely by thermal motions 
(resulting in a temperature r D LA = (1.24 ± 0.06) x 10 4 K), gives 
Ax 2 = 62.0 and is thus ruled out at more than 6.5<x. 

The effect of raising the gas temperature from the best-fitting 
low values to the 2a upper limit is most pronounced for the C 11 col- 
umn density and corresponding C abundance. The reasons are two- 
fold. First, the two Cn lines accessible, /11334 and /11036, are both 
strong and close to saturation (see Figure[TJ. For all other species 
transitions are available which are so weak that the column density 
is almost independent of the Doppler parameter, and is thus insen- 
sitive to the DLA temperature. Second, C is the lightest among the 
metals observed, and is thus the one with the largest value of b th 
at a given temperature. Consequently, adopting the 2cr upper limit 
Tola < 4700 K has the effect of reducing [C/Fe] by 0.37 dex rela- 
tive to the value in Table [T] diluting the degree of C enhancement 
deduced for this DLA. On the other hand, [O/Fe] is unchanged (for 
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Figure 3. On the basis of the available metal lines, we show the differential 
X 2 value, Ax 2 ( t dla) = X 2 ( t dla) -*^ in . Ior a ran g e 01 plausible kinetic 
temperatures for the new candidate CEMP DLA (solid curve). We also show 
the lcr and 2<r confidence regions for a^ 2 distribution with 7 degrees of 
freedom (horizontal dashed lines; see text for further details). This suggests 
that the DLA's kinetic temperature is 7Yjla < 4700 K. A purely thermal 
model, indicated by the dotted lines, is ruled out at 6.5<x. 



the reasons just explained: both O and Fe have weak absorption 
lines in our spectrum); thus, even for the upper limit on the tem- 
perature, the [O/Fe] abundance in this DLA is ~ 2a away from the 
values typical of the metal-poor DLA population (Figure[2j- 

The usefulness of the analysis described above is limited by 
the fact that the masses of the elements included in our data set 
span less than a factor of five between C and Fe, so that the cor- 
responding values of for a given T differ by little more than 
a factor of two (recall that b lb cc /n~ 1/2 ). Further progress could 
be made by observing a much lighter element, such as D. How- 
ever, absorption lines from the Lyman series of Di can only be 
resolved from the neighbouring Hi transitions in high order Ly- 
man lines (e.g. lO'Meara et al.l2006l;|Pettini et al.l2008bl) : recording 
these transitions with the required S/N ratio will involve dedicated 
and time consuming observations in the ground-based UV portion 
of the spectrum. 

In conclusion, while the most likely interpretation (i.e. the one 
producing the lowest value of x 2 ) of the absorption lines in the 
DLA considered here implies a C enhancement by a factor of ~ 4 
([C/Fe] = +0.59), the uncertainties introduced by the range of gas 
temperatures allowed by the data lead us to consider the z abs = 
3.067 DLA in front of J1358+6522 as a candidate CEMP DLA. 



4 GAS TEMPERATURE ESTIMATES IN DLAS 

At present, there have been very few measurements of the kinetic 
temperature for high redshift DLAs. This is hardly surprising, since 
the contribution to the total metal-line broadening in a gas with a 
temperature < 10 4 K is nearly undetectable, even for systems with a 
relatively low velocity dispersion (a ~ 5 km s~'). Indeed, the sole 
direct measurement of a DLA's kinetic temperatur e, on the basis 
of me tal-lines alone, was reported only recently by ICarswell et al.l 
( 120121) . who measured a gas temperature of T DLA = (1.2 ± 0.3) x 
10 4 K for the z abs = 2.076 DLA towards Q2206 - 199. 

For the DLA we report herein, the upper limit on the ki- 
netic temperature of T DLA < 4700 K is consistent with the gas 
residing in either a cold (a few xlO 2 K) or a warm (a few xlO 3 
K) neutral medium. In fact, the absence of molecular absorption 
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lines (or species such as C i) might suggest that the DLA's gas is 
redominantly warm, with a temper ature that is likely > 3000 K 
Petitiean. Srianand. & Ledouxl200dh . As a side note, we also com- 
ment that our derived limit for this system is consistent with the 
best-fitting temperature of 4000 K for the CEMP DLA towards 
J0035 - 0918, using the metal absorption lines alone. 

At present, the majority of DLA temperature estimates rely 
on measuring the spin temperature from the Hi A2\ cm line 
and assuming that it is thermalised (e.g. iKanekar et al.l [20091 ; 
ISrianand et al]|2012h . This may not be a valid assumption in the 
metal-poor regime, where there are fewer atomic coolants avail- 
able to lower the temperature of the gas below T DLA < 1000 K; for 
temperatures above ~ 1000 K, collisio ns are less e fficient at ther- 
malising the Hi A21 cm line (see e.g. lLiszlll200lh . resulting in a 
spin temperature that is somewhat less than the kinetic tempera- 
ture. Although this remains to be tested in a system where both the 
spin and kinetic temperature can be derived, the current upper lim- 
its based on line-profile techniques are consistent with the lower 
limits afforded by spin temperatures for [Fe/H]< -2.0. 

Certainly, the simplicity of the absorption lines from such qui- 
escent DLAs highlights them as ideal systems to directly measure 
the kinetic temperature of neutral gas in the metal-poor regime. 
The ability to measure or place limits on the kinetic temperature 
of DLAs is an exciting prospect for future research; one could in 
principle test how the gas temperature varies as a function of metal- 
licity, H i column density, and redshift. 



5 COMPARISON WITH NUCLEOSYNTHETIC YIELDS 

We now compare the best-fitting abundance pattern shown in Fig- 
ure[2]to illustrative models of stellar nucleosynthesis. The ultimate 
aim of this work is to decipher the clues into the nature of the stars 
responsible for the metal enrichment of the most metal-poor DLAs; 
here we look critically at what advances can be made towards this 
goal using the new data presented here. 

At very low metallicities it is likely that the dominant metal 
yields are from massive stars (10-100 M G ) that explode as core- 
collapse supernovae (CCSNe). At present, however, we still lack 
a full physical understanding of the explosion mecha nism that op- 
erates dur ing the collapse of a massive star (see e.g. iMezzacappal 
120051 and lOtt. O'Connor. & Dasguptall20lH for recent reviews on 
this topic). Despite these uncertainties, several groups have devel- 
oped theoretical models that parameterise the unknown physics, 
such as the degree of mixing between the stellar layers and the final 
kinetic energy of the explosion, in order to estimate the resulting 
nucleosynthetic yields. Here we consider two CCSNe model suites 
to compare to the measured abundance pattern of the candidate 
CEMP DLA: (1) the recent calculations for metal-free Populatio n 
III stars in the mass range 10-100 M G by|H eger & Woosle vtd2010h; 
and (2 ) the Population III/II/T calculations by Ichieffi & Limongil 
(2004) for stars with masses in the range 13-35 M Q . 

The recent database of zero-me tallicity nucleosynthesis calcu- 
lations bv lHeger & Wooslevl d201 d) contains a suite of 1440 mod- 
els, with two parameters describing the physics of the explosion 
and one parameter describing the mixing. Although these calcula- 
tions use a fine mass-resolution of ~ 0. 1 M n , the model stars that 



Table 2. Summary of best-fitting model yields 



Model 


Z(Z Q ) 


Mass (Mq) 


[C/Fe] 


[O/Fe] 


Fe Mass (M Q ) 


HW10 

CL04 

LC12 


0.00 
0.05 
0.00 


17!s 
35 
35 


+0.60 
+0.30 
+0.34 


+0.67 
+0.68 
+0.80 


0.017 6 
0.215 
0.144 



differ in mass by just ~ 0. 1 M Q often produce notable differences 
in the final nucleosynthetic yields. We have therefore opted to in- 
tegrate the yields over a Gaussian-like initial mass function (IMF; 
two additional model parameters are used to define the centroid 
[Mimf] and width [<ri MF ] of the IMF). In this way, we hope to ob- 



a HW1 0:lHeger & Wooslevl feOld) ; CL04: IChieffi & LimongiH2004 ; 
LC12: lLimongi & Chieffil i2012h . 

b For the HW10 model, this is the IMF-weighted ejected Fe mass. 



tain an approximate handle on the mass range of the stars responsi- 
ble for the DLA's enrichment. Thus, the resulting parameter space 
contains 196000 combinations of the five adjustable parameters. 

To determine the best-fitting model, w e use the star_fit soft- 
ware provided by iHeger & Wooslevl d201Ch Fl The combination of 
parameters which best fits the DLA element abundances has an 
IMF centred on Mimf = 17 M Q , with a width of logo"iMF/M G = 
0.15 (details are summarised in Table |2j. This model has a rela- 
tively low degree of mixing between the stellar layers (approxi- 
mated by a running boxcar filter with a width equal to 2.5% of the 
total He core mass), and a final kinetic energy of the ejected mate- 
rial of 1.4 x 10 51 ergs -1 . For this IMF- weighted model, the average 
star will eject 0.017 M G of Fe. An Fe yield of just 0.017 M Q is quite 
low for such an explosion energy, this occurs as a result of the low 
degree of mixing between the stellar layers during the explosion. 

The abundance pattern generated with this model is shown in 
the left panel of Figure|4] Overall, there is good agreement between 
calculated and measured abundances, except possibly for Al which 
deviates by ~ 3cr. However, the discrepancy is only 0. 13 dex, which 
is small considering the likely uncertainties in the calculations of 
model yields. 

We next consi der the nucleosynthesis ca lculations of Popula- 
tion III/II/I stars by |Chieffi& Limongil d2004l) . This database con- 
tains a set of six stars of different mass, spanning a range of metal- 
licity (a further six choices of metallicity are available, including: 
Z = 0.0,5 x 10~ 5 , 0.005, 0.05, 0.3, and 1.0 Z G ), and an effectively 
continuous range of 'mass cut', which is a parameterisation for the 
amount of material that escapes the binding energy of the star fol- 
lowing the explosion; this leads to a model with only three free 
parameters. The model that comes closest to reproducing the chem- 
ical composition of our candidate CEMP DLA is for a 35 M Q star 
with a metallicity of Z = 05 Z a ejecting 0.206 M G of 56 Ni (the 
default value considered by Ichieffi & Limo ngi 2004 is 0. 1 M G of 
56 Ni). A summary of the model details is presented in Table [2] In 
order to eject 0.2 M G of 56 Ni, the kinetic energy of the ejecta for 
a 35 M G star would need to be an order of magnitude larger than 
a typical supernova, a nd will thus reside on the hypernova branch 
(see e.g. Figure 1 from lTominaga. Umeda. & Nomotoll2007|) . 

The abundance pattern for this model is shown in the right 
panel of Figure [4] as the blue stars. Qualitatively, the match to the 
data is poorer than in the previous case, with both C and N appear- 
ing somewhat discrepant. However, we caution against reading too 
much into such a comparison: the abundances of C and N have the 
largest errors among the elements we observed and this model in- 
volves fewer free parameters. Statistically, the two models shown in 



star_fit is available for download from 
http : //homepages . spa . umn . edu/ -alex/znuc/ 



© 201 1 RAS, MNRAS 000,[T]{8] 



A new candidate CEMP DLA 1 




Al 



Fe 



El, 



X 



+ 1.0 






+ 0.5 






0.0 






-0.5 




it 


-1.0 




ft 



-fr Z=0.00 z E 
it Z = 0.05 Z r 



Figur e 4. For the most probable abundance pa ttern (see Figure [2), we o verplot the best-fitting models of Population III (left panel stars; iHeger & Wooslevl 
20 id) and Population II ( right panel blue stars: Ichieffi & L imongi 2004) nucleosynthesis. We also plot the best-fitting Popula tion III model from the recent 



Limongi & Chieffi 



Chieffi & 



121 



yield calculations (right panel orange stars). Statistically, the metal-free IHeger & WoosleM 1201 Ot) model and the Z = 0.05 Z Q 
i model provide an equally good fit to the data, given that the former contains more model parameters. The recent metal-free calcu- 



lations by 
represented by all three models. 



Limongi & Ch ieffi 1 20 1 2), on the other hand, are less favoured. Given the (unknown) model uncertainties, we consider the data to be equally well 



Figure[4]are equally probable, if the (unknown) model uncertainties 
are comparable in the two cases. 

To determine which of the two model suites provide a closer 
representation of the data we use the Bayes evidence ratio (oth- 
erwise known as the Bayes factor). This measure has the advan- 
tage of being able to directly compare two model suites with dif- 
ferent physical prescriptions and a different num ber of free para m- 
eters. Using the Laplace approximation (see e.g. Mackay 2003, p. 
341) and assuming flat pr iors, we measure the e vidence ratio of the 
Ichieffi & Limongi! i2004l) models relative to the lHeger & Wooslevl 
d2010l) models to be £ ~ 1.5, confirming our earlier conclusion that 
it is difficult to discern between the two models. 

In the right panel of Figure [4] we also show the best-fitting 
model from the recent iLimongi & Chieffil J2012h metal-free nu- 
cleosynthesis calculations (orange stars). These new models were 
computed with the latest version of their stellar evolutionary code, 
and now include higher mass stars (with a total mass range of 13M Q 
to 80M o ). The best-fitting model with this suite of calculations is 
for a 35M G metal-free star that ejects 0.139M o of 56 Ni. Although 
the Z = 0.05 Z model gives a somewhat better fit, we consider 
both models to be adequate representations of the data given the 
unknown model uncertainties. 

We have not considered models of rotating low metallic ity 
Population II l lMevnet. Ekstrom, & MaedeJI 20061 ' iHirschil 200% or 
Population III stars l ljoggerst et alj|2010h . These models tend to 
yield a large [N/Fe] abundance, something that we do not observe 
in the abundance pattern of this DLA. We also refrain from comput- 
ing full-scale galactic che mical evolution models. Such an exercise 
was performed recently by I Salvador! & Ferraral | |2012|) to infer the 
enrichment history of CEMP DLAs. These authors conclude that 
such DLAs are enriched by Population III stars, but may receive 
an additional contribution from AGB stars - which would increase 
both the C and N abundances (see e.g. Table 6 from lCristallo et all 
1201 lh . Given the relatively low N abundance observed here, we 
suggest that the metal contribution from AGB stars for this system 
is minimal. 

Thus, on the basis of the available data, it is not possible to 
distinguish statistically between the two possibilities considered 
here, enrichment solely by metal-free stars, or by stars of moder- 
ate metallicity, with Z ~ 1/2OZ (and of course it would be even 
more difficult to recognise a mixed scenario, involving more than 
one stellar generation). Looking ahead, there are two foreseeable 



avenues that may yield tighter constraints on the nature of the stars 
that enrich the most metal-poor DLAs: 

(1) Improvements in the modelling of massive star nucleosyn- 
thesis and, in particular, a better description of the physics behind 
the explosion mechanism. Much progress is currently b eing made 
on the relevant theory (see e.g. iBelczvnski et "all 1 20 111) . Further- 
more, future increases in computational power will allow the con- 
struction of large databases of finer-grid nucleosynthesis calcula- 
tions. This may remove the 'discreteness' of the current models. 

(2) Stronger discriminatory power between different mod- 
els, particularly regarding the explosion energ y, lies with the ra- 
tio of two Fe-peak elements (e.g. rCr.Ni.Zn/Fel. lUmeda & Nomotol 
l2002h . In the most metal-poor DLAs, absorption lines from these 
less abundant elements are too weak to be detected with current ob- 
servational facilities, but may be within the grasp of the next gen- 
eration of 30-m class telescopes. This will open up new avenues 
to probe the details of the explosion. Furthermore, the greater light 
gathering power of these future facilities will make it feasible to 
observe the faint background quasars with the spectral resolution 
R > 100000 required to fully resolved the metal absorption lines 
of most DLAs, and thus decouple directly (rather than by indirect 
means as done here) the relative contributions to the line widths 
from thermal and turbulent broadening. 



5.1 Identifying more CEMP DLAs 

At present, the discriminatory power between Population II and 
Population III models rests with the [C/Fe] ratio. We should there- 
fore compile a list of candidate systems, such as the one described 
here, that exhibit the following hallmark signatures: (1) an en- 
hanced [C/Fe] ratio with relatively normal [N/Fe] (i.e. < 0.0). Cur- 
rent models of Population II nucleosynthesis have difficulty in re- 
producing these signatures. (2) an [O/Fe] ratio that is enhanced be- 
yond t he typical value observed for DLAs (~ +0.4: ICooke et al.l 
l2011bl) : enhanced [C/Fe] is more likely to yield enhanced [O/Fe] 
as well (assuming the Popu lation III nucleosynthesis calculations 
by IHeger & Wooslevl l20ld) . Furthermore, there are several weak 
O i lines that are insensitive to the details of the cloud model and 
in particular to the balance between thermal and turbulent broaden- 
ing. The more easily measured O/Fe ratio may therefore prove to 
be a reliable indicator of further examples of CEMP DLAs. 
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6 SUMMARY AND CONCLUSIONS 

We have reported new observations, at a higher spectral resolution 
than any published previously, of the very metal-poor damped Lya 
system at z abs = 3.067 towards the QSO J1358+6522. Our Keck 
HIRES spectrum covers 21 absorption lines from ion stages of H, 
C, N, O, Al, Si, S, and Fe which are dominant in neutral gas. From 
the analysis of these data we draw the following conclusions: 

(i) The metallicity of the DLA, as measured by the Fe/H ratio, 
is ~ 1/700 of solar ([Fe/H] = -2.84). Both C and O appear to be 
enhanced by a factor of ~ 4 relative to Fe. 

(ii) The gas is 'cold', in the sense that the profiles of the absorption 
lines are best reproduced by a model that has a low value of the 
Doppler parameter, b tot = 3.0 km s -1 , which furthermore appears to 
be dominated by turbulent, rather than thermal, motions. 

(iii) We derived a 2cr upper limit on the DLA's kinetic tempera- 
ture of Tola < 4700 K. The thermal broadening that corresponds 
to this upper limit lowers the [C/Fe] ratio by 0.37 dex, but does not 
affect the [O/Fe] which is deduced from weak, unsaturated absorp- 
tion lines. 

(iv) Finally, we compared the abundance pattern for this candidate 
CEMP DLA to model nucleosynthesis calculations for Population 
I, II, and III stars. Given the handful of elements available and the 
present limitations of the models, we are unable to formally distin- 
guish between models of Population II and Population III enrich- 
ment. Upcoming facilities will permit the abundances for several 
Fe-peak elements to be measured, providing additional diagnostics 
between these two possibilities. 

CEMP DLAs might be more common than currently appre- 
ciated; the available C n lines are almost always saturated, making 
it difficult to measure accurately the C abundance unless one can 
determine a secure cloud model from a careful analysis of the line 
profiles. Given the potential of CEMP DLAs for probing nucle- 
osynthesis by Population III stars, a sample of the most promis- 
ing candidates should now be compiled. With the next generation 
of 30-metre class telescopes, it should become feasible to record 
absorption lines from the intrinsically less abundant Fe-peak ele- 
ments, and thereby explore in more detail than is possible at present 
the properties of the core-collapse supernova explosions of metal- 
free stars. 



ACKNOWLEDGEMENTS 

We are grateful to the staff astronomers at the Keck Observa- 
tory for their expert assistance with the observations. We also 
thank an anonymous referee who provided some comments on the 
manuscript. Valuable advice and help with various aspects of the 
work described in this paper was provided by Bob Carswell, and 
Sergey Koposov. We thank the Hawaiian people for the opportu- 
nity to observe from Mauna Kea; without their hospitality, this 
work would not have been possible. RC is supported by a Re- 
search Fellowship at Peterhouse College, Cambridge. MTM thanks 
the Australian Research Council for a QEII Research Fellowship 
(DP0877998). 



Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, ARA&A, 47, 481 
Becker G. D.. Sargent W. L. W.. Rauch M. Calverley A. P., 201 1, arXiv. 

|arXiv:1101.4399l 
Beers T. C, Christlieb N„ 2005, ARA&A, 43, 531 

Belczynski K., Wiktorowicz G., Fryer C, Holz D., Kalogera V, 2011, 
larXiv:1110.1635l 

Carswell, R. F., Jorgenson, R. A., Wolfe, A. M., & Murphy, M. T. 2011, 

MNRAS, 411,2319 
Carswell R. F, Becker G. D., Jorgenson R. A., Murphy M. T., Wolfe A. M., 

2012, MNRAS, submitted 
Chieffi A., Limongi M., 2004, ApJ, 608, 405 

Cooke R., Pettini M., Steidel C. C, Rudie G. C, Jorgenson R. A., 2011a, 
MNRAS, 412, 1047 

Cooke R., Pettini M., Steidel C. C, Rudie G. C, Nissen P. E.. 201 lb, MN- 
RAS, 417, 1534 

Cooke R., 201 lc, PhD Thesis, University of Cambridge. 

Cristallo S., et al., 2011, ApJS, 197, 17 

Ellison, S. L., Kanekar, N., Prochaska, J. X., Momjian, E., & Worseck, G. 

2012, MNRAS, in press <arXiv; 1205. 12001 . 
Erni P., Richter P., Ledoux C, Petitjean P., 2006, A&A, 451, 19 
Heger A., Woosley S. E., 2010, ApJ, 724, 341 
Hirschi R., 2007, A&A, 461, 571 

Joggerst C. C, Almgren A., Bell J., Heger A., Whalen D., Woosley S. E., 
2010, ApJ, 709, 11 

Jorgenson, R. A., Wolfe, A. M., Prochaska, J. X., & Carswell, R. F. 2009, 
ApJ, 704, 247 

Kanekar N„ Smette A., Briggs F. H., Chengalur J. N, 2009, ApJ, 705, L40 
Kobayashi C, Tominaga N, Nomoto K., 201 1, ApJ, 730, L14 
Lampton M., Margon B., Bowyer S., 1976, ApJ, 208, 177 
Limongi M., Chieffi A., 2012, arXiv. larXiv: 1202.458 II 
Liszt H., 2001, A&A, 371, 698 

Lucatello S., Tsangarides S., Beers T. C, Carretta E., Gratton R. G, Ryan 
S. G., 2005, ApJ, 625, 825 

Mackay D. J. C, 2003, Information Theory, Inference, and Learning Algo- 
rithms. Cambridge Univ. Press, Cambridge 

Meynet G., Ekstrom S., Maeder A., 2006, A&A, 447, 623 

Mezzacappa A., 2005, ARNPS, 55, 467 

O'Meara J. M., Tytler D., Kirkman D., Suzuki N., Prochaska J. X., Lubin 

D., Wolfe A. M., 2001, ApJ, 552, 718 
O'Meara, J. M., Buries, S., Prochaska, J. X., et al. 2006, ApJ, 649, L61 
Ott C. D., O'Connor E. P., Dasgupta B., 2011, larXiv:1111.6282l 
Penprase B. E., Prochaska J. X., Sargent W. L. W., Toro-Martinez I., Beeler 

D. J., 2010, ApJ, 721, 1 
Petitjean P., Srianand R., Ledoux C, 2000, A&A, 364, L26 
Pettini M., 2004, in Esteban C, Herrero A., Garcia-Lopez R., Sanchez F., 

eds, Cosmochemistry: The Melting Pot of Elements. Cambridge Univ. 

Press, Cambridge, p. 257 
Pettini M., Bowen D. V, 2001, ApJ, 560, 41 

Pettini M., Zych B. J., Steidel C. C, Chaffee F. H., 2008a, MNRAS, 385, 
2011 

Pettini, M., Zych, B. J., Murphy. M. T., Lewis, A., & Steidel, C. C. 2008b, 

MNRAS , 391, 1499 
Ryan S. G., Aoki W., Norris J. E., Beers T. C, 2005, ApJ, 635, 349 
Salvadori S., Ferrara A., 2012, MNRAS, 421, L29 

Srianand R., Gupta N, Petitjean P., Noterdaeme P., Ledoux C, Salter C. J., 

Saikia D. J.. 2012, MNRAS, 421, 651 
Tominaga N, Umeda H., Nomoto K., 2007, ApJ, 660, 516 
Tumlinson J., et al., 2010, ApJ, 718, L156 
Umeda H., Nomoto K., 2002, ApJ, 565, 385 
Umeda H., Nomoto K., 2003, Nature, 422, 871 
Vogt S. S., et al., 1994, SPIE, 2198, 362 

Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 2005, ARA&A, 43, 861 



REFERENCES 



Aoki W., Beers T. C, Christlieb N., Norris J. E., 
S., 2007, ApJ, 655, 492 



Ryan S. G., Tsangarides 



© 201 1 RAS, MNRAS 000.[T 



